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Abstract 

If the Standard Model (SM) of particle interactions is extended to include a second scalar doublet 
[H^ , [H^ + iA^)/\/^], which is odd under an unbroken Z2 discrete symmetry, it may be called the 
dark scalar doublet, because its lightest neutral member, say H^, is one posssible component for 
the dark matter of the Universe. We discuss the general phenomenology of the four particles of 
this doublet, without assuming that is the dominant source of dark matter. We also consider 
the impact of this dark scalar doublet on the phenomenology of the SM Higgs boson h. 
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I. INTRODUCTION 



The canonical Standard Model (SM) of quarks and leptons requires only one Higgs doublet 
{(f)'^ , (f)^) , such that (fp acquires a nonzero vacuum expectation value {(jP) = v = 174 GeV, 
from which all particles (except neutrinos) obtain mass. Suppose a second scalar doublet 
[H^, + iA°)/-\/2] is added [l], which is odd under an unbroken Z2 discrete symmetry, 
then there are at least two interesting outcomes. 

(I) This scalar doublet also has a conserved additive quantum number, so that {H^ + 
iA^)/^/2 is a mass eigenstate, i.e. tuho = m^o. Such is in fact the case of supersymmetry, 
where this is nothing but the complex conjugate of a scalar lepton doublet. More generally, 
it may belong to a new class of particles [2, yj with a separately conserved additive quantum 
number unrelated to lepton number. 

(II) There is no additional symmetry beyond the postulated Z2, in which case rrifjo 7^ m^o. 
As pointed out [l| already in 1977, either or must then be stable. This idea has been 
revived recently, first in the context of radiative neutrino masses [4], then as a means of 
understanding precision electroweak data [5]. In any case, this specific potential new source 

is gaining more attention [?, 8, 9, 



101, 
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12|. It will also 



of dark matter of the Universe 
have an important bearing on interpreting forthcoming data from the LHC (Large Hadron 
Collider) at CERN. 

In this paper we will analyze in detail the two-scalar-doublet structure of scenario (II), 
which we call DSDM (Dark Scalar Doublet Model). Whereas there are specific constraints 



from identifying if ° as the sole source of dark matter 



0,3, 



we will take the viewpoint that 



there may also be other sources of dark matter such as that coming from supersymmetry. For 
example, the first such explicit model with two (and possibly three) different dark-matter 
particles has already been proposed flS]. Thus our analysis will allow a wider range of 
parameter values, with the aim of extracting the possible signals of this dark scalar doublet 
at the LHC. This doublet has been called the "inert Higgs doublet" , but it is neither 
inert (since it has gauge and scalar interactions) nor a contributor to the Higgs mechanism 
(since it has no vacuum expectation value), hence we prefer to call it the dark scalar doublet. 

Consider the scalar potential of $1 = (0+,0°) and $2 = + iA'^)/\^], with $1 

even and $2 odd under an exactly conserved Z2: 

V = /i?$I$i+/i2$t$2 + lAi($t$i)' + ^A2($^$2)' + A3($I$l)($|$2) 



+ A4($I$2)(*5*l) + ^A5($I$2)' + ^X*M<^lf. (1) 

All the above parameters are necessarily real except A5, but the exact Z2 symmetry (which 
forbids the bilinear term $|$2 allows us to rotate the relative phase between $1 and $2 
to make A5 real as well. As the standard SU(2) x f/(l) gauge symmetry is spontaneously 
broken by (0°) = v = 174 GeV, the masses of the resulting scalar particles are given by |j] 

m'^{h) = 2Al^;^ (2) 

m\H^) = fil + Xsv', (3) 

m^H"^) = ^1 + {X, + X, + X,)v\ (4) 

m\A') = + (A3 + A4 - X,y. (5) 

The lone Higgs boson of the Standard Model is of course h, whereas H^, H^, and are 
the components of the dark scalar doublet which interact with h and among themselves as 
follows: 



Vint — -X2 



H+H- + ^iHy + ^iAy 
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'H+H~ + ^{Hy + ^{Ay 



+ ^(A4 + A5) (vh + i/.^) [Hy + i(A4 - A5) (vh + {Ay. (6) 

Being electroweak doublets, they also have gauge interactions, but they do not interact 
directly with quarks or leptons. 

Assuming the mass hierarchy mH± > ttlao > 171^0 , the stable if" then appears as missing 
energy in the decays of A^ and H^. Since there is no term linear in A^ or in Eq. (6), 
the decay of A^ must occur through the gauge interaction 

^ Z.(if°9M°- (7) 

2 cos 6w 

Hence the dominant decay of A'^ is into H'^ff, where / = lepton or quark, and is missing 
energy. 

As for if^, its gauge interactions with if" and ^4° are given by 

^W-{H^d>'H+ - + |w7(A°9^/f + - H+d^'A'^) + h.c. (8) 

Hence the dominant decays of H"^ are into W^H^ and ly^A*^, with the subsequent decay 
of A° into Z^H^ in the latter. 

The only other trilinear gauge interactions of the dark scalars are those of iJ^ with Z 
and 7. There are of course also the quadrilinear terms required by gauge invariance. 



II. LEP CONSTRAINTS 



A. Z-boson decay 

If mffo + m^o < mz, then Z can decay into H^A^ with the subsequent decay of A'^ to 
f f where / is a quark or lepton. This would appear as a rare Z decay event with // + 
missing energy in the final state. The high precision LEP-I experiments have accumulated 
a sample of about 1.7 x 10^ on-shell Z boson decays via the qq and t^^~ channels jl^ . 
ALEPH preformed a partial analysis of an integrated 79 pb~^ to search for the rare decay 
Z —>■ ff + uu and found no event above the SM prediction [isl]. To estimate how the whole 
LEP-I data sample can constrain the DSDM, we calculate the SM partial decay width of 
Z ^ f f + uu {f = q, €) and obtain its decay branching ratio 

^^,ry rr , (4.269gg-,p + 0.779,+,-,p) X 10"^ GeV 

BK [Z J J + uu^ — 



2.494totai GeV 
^ 2.02 X 10-^. 

We thus expect about (1.7 x 10'') x (2.0 x 10"'') ~ 3.4 background events from the SM. 
Assuming that only 4 events show up as // + missing energy from the whole LEP-I data 
sample, we deduce according to Possion statistics that only about 5 signal events are allowed 
at 95% C.L. in the ff + H^H° channel, i.e. 

^ ^ < ~ 3.7 X 10"^ 

To ~ 1.7 X 10^ X 0.8 

where the SM branching ratio BR[z — ^ //^ ~ 0.8 (for / = q,tj has been used. This 

implies that the phase space for Z —>■ H^A^ is very small (with a mass gap of a few MeV), 



so for all practical purposes, the bound rriHO + itiao > mz should be respected HI]. The 



resulting lower limits on rrifjo for various Aitlaoho = ^1^0 — tjiho are shown in Table [11 

TABLE I: Lower limits on mjjo from the LEP Z-decay constraint for various Amy^oji/o. 
Am^oi^o (GeV) 5 10 15 20 30 



niHO (GeV) 43 41 38 36 31 



Since A^ — > + vv is also possible, this process may also contribute to the invisible 
width of the Z if m//o -|- m^o < niz- However, the corresponding experimental uncertainty 
allows for a mass gap of a few GeV, which is not competitive with the constraint discussed 
above. 
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B. Neutralino search at LEP 



While there is no existing experimental search that can be directly translated into a 
search for the dark scalars, the search for charginos and neutralinos isp^ry close as they 
may exhibit the same collider signatures. The OPAL Collaboration [16] has performed 
a search for neutralino pair production via the process of e^e^ X2X1 the hadronic 
decay mode X2 ^ Xi^* XiQQ where Xi 2 denotes the lightest (next lightest) neutralino, 
respectively. No significant excess above the SM background leads to an upper limit on the 
cross section of X2X1 associated production: cr~o~o ~ 0.1 pb — 0.5 pb for m~o ~ 50 GeV and 



60 — 80 GeV (see the right figure of Fig. 10 in Ref. [3]). The L3 Collaboration 
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also performed a similar search at y/s = 189 GeV, which gives the upper limit on the cross 
section as o"~o~o ~ 0.1 pb — 2pb for m~o ~ 50 GeV and m^o ~ 60 — 80 GeV (see Fig. 3a in 
Ref. |l7|). 

The translation of the above limits on the cross section from the neutralino to the dark 
scalar is not straightforward even though they share the same topologies. The kinematics 
and spin correlations of X1X2 A^H^ pair productions are totally different. For example, 
the scalar pair has to be in the p-wave as it is produced via a virtual Z-boson, therefore. 



each scalar has large transverse momentum ISj. Furthermore, there is no spin correlation 



between the two scalars. Those differences lead to different acceptances of the kinematics 
cuts used by the OPAL and L3 Collaborations to extract the signal from the background. 
One thus cannot simply apply the OPAL and L3 constraints to the DSDM. In order to get the 
realistic limits, one has to go through the whole analysis using the correct kinematics and spin 
correlations. For a rough estimation, we now simply compare the prediction of the DSDM 
to the OPAL and L3 limits. The production cross section of e^e^ —>■ A^H^ in the DSDM at 
v/i = 208 GeV is about 0.28 pb - 0.20 pb for m^o = 50 GeV and m^o = 60 GeV - 80 GeV. 
It is reasonably consistent with both OPAL and L3 constraints. Hence, we will use these 
masses in the following monte carlo study. 

C. Higgs direct search 

If the Z2-odd scalars are lighter than the SM Higgs boson, they may have a great impact 
on the direct search of the latter because h H^H^ and h — > A^A^ may become the 
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FIG. 1: (left) The 95% C.L. upper limit on BR{h 



invisible), adapted from [19]; (right) The 



20( 1 ■ The label H in these figures refers to the 

2 



20| for 



95% C.L. upper limit on BR{h — > 66), adapted from 

SM Higgs boson h, and in the right figure is defined as = {dHZZ / gfi^zz) (^^^ 
details) . 



dominant decay channels. In particular, the h H^H^ channel is invisible. LEP II has 
performed a direct search of the SM Higgs boson in its invisible decay mode but found no 
event, thus obtaining an upper bound on this decay branching ratio, shown on the left in Fig. 
[H Below we will use this information to find the allowed parameter space of the DSDM. 
As for the h — » A^A^ decay mode, it may give rise to four charged leptons plus missing 
energy, resulting from A^ — > H^Z* H^i^i^. At LEP the SM Higgs boson is produced 
in association with an on-shell Z-boson, i.e. e~^e^ Zh, which decays into two charged 
leptons or two jets. Thus, if the SM Higgs boson decays into two A^ bosons, there could 
be a collider signature of six charged leptons plus missing energy, or four charged leptons 
plus two jets plus missing energy. A new scan over the accumulated data sample at LEP to 
search for such a signature may prove to be interesting. 

LEP II searched the SM Higgs boson mainly via the process e~^e~ Z{-^ £~^£~)h{—>- bb) 
as the decay h bb dominates for a light Higgs boson. The negative search at LEP II gives 
rise to the upper bound on the decay branching ratio oi h ^ bb, shown in the right in Fig. 
[11 and sets a lower limit of the SM Higgs boson mass of nih > 114.4 GeV. However, this 
limit can be relaxed in the DSDM because the invisible decay of the Higgs boson may be 
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dominant so that the decay branching ratio into bb is highly suppressed ^. 

The LEP II bounds on h invisible and h ^ bb compete with each other. The former 
requires smaller invisible decay branching ratio while the latter requires larger. Consider 
now the possibility of the DSDM satisfying both limits. 

The partial decay width of hsM bb is given by 



h 



while the partial decay width of hsM SS {S denoting a member of the dark scalar 
doublet) is given by 



Y(h^ SS) = 5s- A^. 

IQrcmh \ 



Arrig 
mi 



(9) 



where Sho/ao = 1 and 6h+ = 2, and the coupling \s is given by 

A//0 = A3 + A4 + A5, Aao = A3 + A4 — A5, Xh+ = A3. (10) 
From Eqs. ([MSD, the above couplings can be further simplified as 

■^Ho = 5 , Ayio = , Ah+ = 5 . [i-i-) 

V V V 

Thus, the partial decay width of the SM Higgs boson, given in Eq. [9l only depends on two 
parameters, ms and fi2, for a given rrih. We note that /i2 is not a free parameter because of 
the requirement of vacuum stability, i.e. 



Ai,2>0; As, A3 + A4- IA5I > -VA1A2. (12) 
Hence we obtain the following bound on /i2, 

fil < mjjo + ]j^mhv, (13) 

which reveals the fact that the self-interaction of the dark scalar doublet, i.e. A2, has a 
nontrivial impact on the SM sector. As to be shown later, the stronger the dark sector 
interacts, the larger is the allowed parameter space of the DSDM for the SM Higgs boson. 

Assuming the charged Higgs boson mass is much larger than rrih and m^o— m^^o = 10 GeV, 
we now scan over the parameter space of the DSDM. In summary, the following conditions 
need to be satisfied in order to relax the LEP limit rrih > 114.4 GeV: 



Similar studies in various models have been carried out in Refs. 



21 



22 



23 



24 



25| 
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LEP II bound, BRmax {inv), on h ^ invisible 



BR (inv) = — 7 ) ^ BRmnr (inv) 



r [hsM ^bb)+T {hsM ^ nm^) + r {hsM ^ a^a^) 

LEP II bound, BRmax {bb) , on h ^ bb 
BR (bb) = —, ^ < BRmax (bb) 



(14) 



ri 


^ bb) 


1 


ri 


[hsM bb) 


1 + r {hsM ^ Hm^) + T {hsM ^ AMO) 



(15) 

• Z-decay constraint 

< mi:/o + m^o. (16) 

• Vacuum stability bound given in Eq. [131 

We note that the vacuum stability bound is highly correlated with the bound on the decay 
of h bb. It can be understood as follows. Due to the Z-pole constraint nifjo should be 
larger than 41 GeV for AniAOHO = 10 GeV, therefore, the SM-like Higgs boson has to be 
heavier than 82 GeV in order to open the h H^H^ mode. Bearing in mind the fact that 
BR^ax (bb) < 0.1 for a light SM-like Higgs boson, e.g. nih ~ 100 GeV, the LEP II bound 
on h ^ bb now becomes 

< BK,^ (») { ^^^s^) ' (17) 



where = y 1 — 4m^/m| and ~ \ ^ ^ ^m^o/w^l- As ~ for a light Higgs boson, 
/i2 has to be very large, even larger than mno, in order to satisfy the above inequality. On 
the other hand, /i2 is also bounded from above due to vacuum stability. It thus leads to 
an interesting connection between the vacuum stability bound and the LEP II bound on 
h bb. 

The allowed region of A2, which satisfies the constraints on both BRmax {inv) and 
BRmax is shown in Fig. [2]^a). For the sake of illustration, we plot the allowed pa- 
rameter space of {niHo, 1^2) for nih = 109 GeV and 112 GeV in Fig. [2](b) and (c). In the 
DSDM the lower limits of nih are shown in Table HTl for various Amyioj/o. 
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FIG. 2: (a) Allowed parameter space (gray region) of {mh,\2)] (t))(c) allowed parameter space of 
{rrifjo, 1-12) after imposing the LEP II constraints on the decay branching ratio of /i ^ invisible and 
h — > bb. The region above the dashed line is excluded by the vacuum stability requirement. From 
top to bottom: A2 = 1.0, 0.5,0.3, O.land A2 0. 

TABLE II: Lower limits on the mass of the SM Higgs boson derived from LEP II Higgs boson 
search for various Am^oj^^o. 

Am^oj^o (GeV) 5 10 15 20 30 
rrih (GeV) 105 106 109 110 110 

III. IMPACT ON THE SM HIGGS SEARCH AT COLLIDERS 



At the LHC, the SM Higgs boson is mainly produced through gluon-gluon fusion induced 
by a heavy (top) quark loop. Once produced, it will decay into a fermion pair or vector 
boson pair. The strategy of searching for the Higgs boson depends on how it decays and 
how large the decay branching ratio is. We present selected decay branching ratios of the 
SM Higgs boson as a function of rrih in Fig. [3]^a), as well as the total decay width {Tfot) 
in units of GeV. If the SM Higgs boson is lighter than 130 GeV, it decays mainly into a 
bottom quark pair (bb). Unfortunately, it is very difficult to search for the Higgs boson 
in this mode due to the extremely large Quantum Chromodynamics (QCD) background at 
the LHC. However, the /i ^ 77 mode can be used to detect a Higgs boson with the mass 



below 150 GeV [26|, 1271] even though the decay branching ratio of this mode is quite small, 
~ O(10~^). If the Higgs boson mass {rrih) is in the region of 130 GeV to 2'mz, the h ZZ* 
mode is very useful because of its clean collider signature of four isolated charged leptons. 
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FIG. 3: (a) Selected Higgs boson decay branching ratios as a function of in the SM; (b) 
selected SM Higgs boson decaying ratios as a function of m/j in the DSDM. Here we have chosen 
miio = 50 GeV, Am^o//o = 10 GeV, mfj± = 170 GeV and ^2 = 20 GeV. The vertical axis is units 
of GeV for the total decay width. 

The h WW^*'^ mode is also important in this mass region because of its large decay 
branching ratio [28]. if mh > 2mz, the decay mode h ZZ —>■ is considered 

as the "gold-plated" mode which is the most reliable way to detect the Higgs boson up to 
rrih ~ 600 GeV because the backgrounds are known rather precisely and the two on-shell 
Z bosons could be reconstructed experimentally. For nih > 600 GeV, one can detect the 
h ZZ i^i^uu decay channel in which the signal appears as a Jacobian peak in the 



missing transverse energy spectrum 29|] . In this section we show a few representive cases to 
illustrate how the new dark scalars affect the Higgs search at the LHC. Our choice of masses 
for all the particles is not necessarily constrained by precision electroweak data [sl, [sol. The 
numerical calculation is carried out with the help of the HDECAY program |3l|, after the 
implementation of the new decay channels in the DSDM. 

The new decay channels of the SM Higgs boson will increase its decay width but decrease 
the decay branching ratios of the usual decay modes, such as bb, WW, ZZ and 77, etc. In 
Fig. [3](b) we present the relevant decay branching ratios of the SM Higgs boson as a function 
of its mass. Clearly, the usual decay modes are highly suppressed once the H^H^ and A^A^ 
decay modes open. The invisible mode {h H^H^) dominates in the light and intermediate 
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FIG. 4: (a) Total decay width of the SM Higgs boson as a function of rrih; (b) ratio of the total 
decay width of the SM Higgs boson in the DSDM and Higgs boson in the SM; (c) decay branching 
ratio of the invisible decay mode of the SM Higgs boson in the DSDM. For comparison, we choose 
mjio = 40 (50, 60) GeV, Am^o/^o = 10 GeV, mjy± = 170 GeV, and ^2 = 20 GeV. 

mass region of rrih, i.e. 100 < < 160 GeV. Due to the 10 GeV mass gap between 
and A^, the contribution of the H^H^ mode is much larger than that of the A^A^ mode, as 
shown by the black solid and dashed curves in Fig. [3](b). It implies that the SM Higgs boson 
decay depends mainly on 111^0 but not m^o. We also note that the 77 mode is suppressed 
so much that it will be very challenging to use this mode to detect the SM Higgs boson in 
the mass region of mh ~ 100 — 150 GeV. 

To examine the dependence of the SM Higgs decay upon irtijo , we choose sample points of 
mj^o being 40 (50, 60) GeV and present the corresponding Higgs boson total decay width, 
the ratio of total decay width in the DSDM to the one in the SM {T /Tsm) and the decay 
branching ratio of the invisible decay mode {Br {invisible)) in Fig. |l](a) (b) and (c), respec- 
tively. Clearly, the total decay width increases by a big amount once the new decay channel 
opens. For example, as shown in Fig. IH^b), the new invisible decay mode can enlarge the 
total width by a factor of three for m^^o = 60 GeV and rrih = 130 GeV and by a factor of two 
for rriHO = 40 GeV and nih = 100 GeV. The decay branching ratio of the invisible mode is 
about 50 — 65% in the mass region rrih ~ 100 — 150 GeV for tuho ~ 40 — 60 GeV. Needless to 
say, the WW mode dominates over the invisible mode once the SM Higgs boson is heavier 
than 160 GeV, beyond which the SM Higgs boson in the DSDM is essentially the same as 
that of the SM. 
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FIG. 5: (a) Ratio of the decay branching ratios of hsM 77 in the DSDM and the SM; (b) 
suppression factor of the usual decay modes compared to the SM for rriffo = 40, 50, 60 GeV. Here 
we have chosen Arrij^ofjo = 10 GeV and fi2 = 20 GeV. 



In the DSDM the decay mode of /i — 77 receives an additional contribution from the 



charged Higgs boson loop. The partial decay width of 77 is given by 



32 



33] 



12872' 



NcQ}9hffAi/2 (Tf) + Qhww-Ai {tw) + ^"r^ A {.rH± 



For the numerical 



where the amplitude Ai and coupling Qhff/hww are given in Refs. [32l, 
calculation we also implemented the correction from the charged Higgs boson loop in the 
HDECAY program. The impact of the invisible decay on the other regular decay modes 
are shown in Fig. [51 We note that the loop correction of the charged Higgs boson slightly 
suppresses the partial decay width of h ■yy by about 10%, as shown in Fig. M^a). Such a 
suppression makes it more difficult to search for the SM Higgs boson via the di-photon mode. 
Furthermore, as an overall effect, all the usual decay modes of the SM Higgs boson, e.g. 
bb/WW/ZZ/'j'j, are highly suppressed once the new decay mode opens. In the intermediate 
mass regime, 130 < rrih < 160 GeV, the suppression factor can be larger than 0.5 which 
definitely makes the search of the SM Higgs boson more challenging. 

On the other hand, one might search for the Higgs boson through its invisible decay 
mode. For example, the weak-boson fusion (WBF) process, qq —>■ q'q'VV q'q'h with the 
subsequential decay to undetectable particles, i.e. h H^H^ /A^A^. The signal is thus 
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FIG. 6: Allowed parameter space of {mjjo,fj,2) to reach the 5o" discovery of the invisible Higgs 
decay via the WBF process: (a) rrih = 110 GeV, (b) mh = 150 GeV. The region above the dashed- 
lines is excluded by the vacuum stability requirement. From top to bottom: A2 = 1.0, 0.5, 0.3, 0.1 
and A2 — > 0. 

characterized by two quark jets, which typically stay in the foward and backward regions 
of the detector and are widely separated in pseudorapidity, and also by a large missing 
transverse momentum {^t), due to the Higgs invisible decay products. It has been shown 
in Ref. that the WBF process can be used to detect the invisible decay Higgs boson at 
5a level, when statistical errors only are considered, if the invisible decay branching ratio 
satisfies the following condition 

as X Br {h invisible) x C 



> 5, 



(18) 



he signal (back- 
s', aB = 167fb 



y/as X C 

where £ denotes the integrated luminosity and asias) the cross section of 
ground) respectively. After imposing the kinematical cuts suggested in Ref. 
and as = 99.4 (99.7, 94.3, 89.2) fb for ruh = 110 (120, 130, 150) GeV. From the above re- 
lation, we obtain the allowed parameter region of {771^0,^2) to reach the So" discovery as 
shown in Fig. [6l Most of the parameter space can be covered for C = 100 fb~^, c.f. the 
yellow region. 

Unfortunately, finding or constraining an invisibly decaying Higgs boson signal in jj + 
event is essentially a counting experiment since a resonance in the invariant mass distribution 
of the Higgs decay products cannot be extracted. It also depends on two or more independent 
model parameters, thus one can only probe or constrain their combination. Next we are 
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going to examine the dark scalar pair production at the LHC, which is complementary to 
the WBF process. 



IV. COLLIDER PHENOMENOLOGY OF THE DARK SCALAR DOUBLET 
A. Production of dark scalar pair at the LHC 

Since all new scalars carry a Z2-odd quantum number, they cannot be produced singly at 
the collider. Furthermore, these new scalars only couple to the Higgs boson and electroweak 
gauge bosons of the Standard Model. Hence they should be produced via processes of the 
Drell-Yan type at the LHC as listed below: 

qq' -^W* ^ A'^H^, (19) 

qq' ^W* ^ H^H^, (20) 

qq^ Z* ^ A^H\ (21) 

qq Z(7) ^ (22) 

In the DSDM, the vertex oi V — (j) — (j) is a pure gauge coupling, where V denotes an 
SM gauge boson while is a Z2-odd scalar. Therefore, the total cross section of each 

production process depends only on the masses of the Z2-odd scalars in the final state. In 
this work the matrix elements of both signal and background processes are calculated using 



MADGRAPH 
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361] while the widths of the dark scalars are calculated in CALCHEP 



37l | with the model file generated by LANHEP [38|]. Agreement of both programs has been 
checked at the level of dark scalar pair production. In Fig. [7] we present the contours of 
the cross sections (in fb units) of the A^H^ and A^H^ associated production at the LHC 
in the plane of the scalars' masses. The red symbol "X" denotes one of our benchmark 
pionts, [mHO, ^ao, fnH±] = [50, 60, 175] GeV, used in the collider study. We note that the 
production cross section of the dark scalar pair at the LHC is sizable, e.g. at the benckmark 
point shown in Fig. [7] the cross section of the A^H^ pair production is about 3pb while that 
of the A^H"^ pair production is about 200 fb. The latter is highly suppressed owing to the 
large mass of H^. Due to the same reason the cross section of the H^H pair production 
is quite small in the relevant parameter space (around the benchmark piont); it is therefore 
not shown in this study. 



14 



200 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I hJ.!, I I I I I I I I I I I 

(a) ■ 



50 




150 - 



100 - 



2e+04 

1e+04 

5e+03 

4e+03 

3e+03 

2e+03 

1e+03 

900 

800 



700 -- 

600 ■■ 

500 ■■■■ 
400 
300 

200 ■■ 

100 - 



20 



40 



60 



80 



100 



100 



90 r 



> 
a> 
CD 



60 ^. 



50 



40 




80 r *'•. 



70 



I i^vi I I* J ^^^ I f.i I I I I i\ I I I I I'l. I \i~ 



40 50 60 70 80 90 100 

(GeV) 

FIG. 7: Contours of the cross sections (in fb units) of A^H^ production (a) and of A^H^ production 
(b) at the LHC in the plane of the scalars' masses. The red symbol "X" denotes our benchmark 
point used in the collider study. Note that Fig. (a) also holds for H^H^ production. 
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FIG. 8: Mass pattern that determines the decay pattern of and H~^. 
B. Decay pattern of and 

The decay pattern of and is completely determined by three mass spacings 
Am^o/^o, Arrifj+fjo and Amfj+j^o, as shown in Fig. [HI where ArriH+fjo = mH+ — mffo and 
Arriff+j^o = 171^+ — m^o. Using the constraint from precision electroweak data, it has been 
shown by Ref. jsl that a relationship exists between the first and second spacing in the 
Z2-odd scalar spectrum 

{mH+ - ruAo) {mH+ - m^o) = M^, M = UOtfo GeV, (23) 

which implies that the charged Higgs boson should be heavier than both A^ and so that 
A^ can only decay into the H^Z* mode. This is consistent with our choice. It corresponds 
to the top-left corner in Fig. [HI Of course, the precision electroweak constraint also requires 
the SM Higgs to be very heavy in this case, say nih > 400 GeV, but that does not affect 
the production and decay of the dark scalar doublet particles which are independent of 
rrifi. [Our view is that the masses of all the scalars being considered here are not necessarily 
constrained by precision electroweak data or dark matter relic abundance because there may 
be other contributions. We simply explore a set of reasonable possibilities.] 

Furthermore, in order to avoid the huge QCD (quantum chromodynamics) background 
at any hadron collider, one has to use the leptonic decay mode to tag the signal event. Thus, 
in Fig. [9] we plot the decay branching ratio (BR) of A^ H^i^i^ and H^i'^u as 

a function of Auiaoho aiid Auih+ho respectively. In these figures we have summed over 
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FIG. 9: Decay branching ratio of the — > H^£~^£~ and the H^£~^u where we have summed 

over two lepton flavors (electron and muon). Note that Arrij^+ffo and Arrij^OffO should satisfy the 
electroweak precision test bound given in Eq. [251 

two lepton flavors (electron and muon). We note that BR (A° H'^i+i") is only sensitive 
to ArriAOHo. If the mass spacing is small, say ArriAOH" ~ lOGeV, the BR is about 0.09. 
As the mass spacing increases, the BR decreases because the massive quark modes, e.g. 
ss{cc,bb), increase. For a very large mass spacing, the BR approaches to the BR of an 
on-shell Z-boson. The decay of is a little complicated because it can decay into two 
modes, H'^W^*^ or A'^W^'-*''. For simplicity, we only focus on the decay mode H^W^*^ 
in this study because the A'^W^*'^ mode involves several particles in the final state 

which render the analysis of the collider signals much more complicated. Again, we note 
that the BR of H^l^v decreases as the mass spacing AmH+H° increases because of 

the massive quark modes. The peaks in the figure are due to the mass threshold of the 
VT-boson. Eventually, the BR approaches the SM VT-boson decay branching ratio, ~ 1/9, 
at the large mass spacing regime. 

In the following sections, we present the results of a detailed Monte Carlo simulation 
and determine the discovery potential of the Z2-odd scalars at the LHC We find that the 
A^H^/H^H^/H^H~ pair productions suffer from huge backgrounds such that the dark 
scalar signals are always overwhelmed. This difficulty persists even if we lower the mass of 
H^. We thus focus on the A^H^ associated production in the next section. 
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C. A^H^ production and signature 



The neutral scalar bosons can be produced via the following process 

pp^Z A°H\ A° ^ H°Z^*\ ^ //, (24) 

which gives rise to a collider signature as //+ ^t- Thus, the dominant SM background 
comes from two processes: 

pp^ ZZ/Z-f, with Z ^ i^u, Zj-i^tr, (25) 

ly+iy-, with iy+ ^ £+z/, w-^rv, (26) 

where two charged leptons in the WW background are required to have the same lepton fla- 
vors. Besides these two intrinsic (irreducible) sources, there are also other reducible sources 
of background such as tt pair production and Wt associated production. For example, if 
the bottom quark from top quark decay escapes the detector, they will contribute to the 
mssing transverse momentum and thus mimic the signal events. But one can suppress these 
reducible backgrounds by vetoing any additional jet activity in the central rapidly region. 
Hence, we focus our attention only on the intrinsic background hereafter. 

Below we choose a benchmark point, [mji^o, m^o, mj:^+] = [50, 60, 170] GeV, to illustrate 
how to discover such a signature at the LHC. In order to mimic the detector, we first require 
both leptons in the final state to satisfy the following basic cuts: 

pfp > 15 GeV, |r/^| < 3.0, (27) 

where p^iv^) denotes the transverse momentum (rapidity) of the charged lepton i. Here 
we have assumed a perfect detector that can precisely measure the four-momenta of the 
final-state leptons. The number of signal and background events, after imposing the basic 
cut, is given in the second column of Table IIIII The number of signal events is much less 
than the one of background events, e.g. the ratio of the signal and background event (S/B) 
is about 10"'^ where S and B denotes the number of the signal and background events, 
respectively. Obviously, additional kinematical cuts are needed to enhance the ratio of 
signal to background. For that, we examine the kinematics differences between signal and 
backgroud below and find out the optimal kinematical cuts to entangle the signal out of the 
background. 

18 




20 40 60 80 100 -1 -0.5 0.5 1 -1 -0.5 0.5 1 

m (GeV) Cos9, , coscj), , 

12 'l 2 12 



FIG. 10: (Normalized) kinematic distributions of the A^H^ production for the benchmark point 
[mfjo,my\^o,mf{+] = [50, 60, 170] GeV (black curve), the WW process (red curve) and the ZZ 
process (green curve) after imposing the basic cuts. The blue dashed-line indicates the optimal 
kinematical cuts which can be used to suppress the SM background. 

In Fig. [TD] we present (normalized) kinematics distributions of both the signal event and 
the background event after imposing the basic cuts, where 9e£{4>u, fnu) denotes the open 
angle (azimuthal open angle, invariant mass) of the two charged leptons and the missing 
transverse momentum originated from the two dark scalars in the final state. We note that 
the kinematics of the signal events significantly differs from the one of background events in 
the following aspects: 

• Since the two charged leptons in the signal event come from an off-shell Z-boson decay, 
they populate the small pt region and they prefer to move parallel to each other, see 
Fig. [TOlfa.b.e.f). On the other hand, the charged leptons of the background events 
have large pr and large open angles. 

• Since the two missing scalars in the signal events are massive, one may expect a large 
^T, but unfortunately, the two missing scalars prefer to move back-to-back and their 
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combination actually results in a small I^t distribution, see Fig. fTOT c). The distribution 
of the ZZ background peaks around 40 GeV while the one of the WW background 
peaks in the large I^t region. 

• The invariant mass distribution of the two charged leptons of the signal events peaks 
around ~ |Am^oj|^o as they come from the decay of an off-shell Z-boson. The ZZ 
background distribution exhibits two peaks: one is around mz (from the on-shell 
Z-boson decay) and the other is around zero (from the virtual photon decay). The 
intermediate region between the two peaks is due to the interference effects between 
the ZZ and Z7 processes. The distribution of the WW background peaks in the large 
invariant mass regime. 

Taking advantage of the kinematic differences between the signal and the background, one 
can impose several optimal kinematic cuts, indicated by the blue dashed-lines in Fig. [TOl to 
extract the signal from the background. 

Based on the topological difference between the signal and background events, we impose 
the following optimal cuts to futher suppress the background, 

< 40 GeV, $T < 60 GeV, 

cos^«>0.9, cos (f)ie> 0.9. (28) 

At last, we require that the invariant mass of the two charged leptons satisfies the following 
mass window cut: 

< m« < 10 GeV. (29) 

In Table IIIIl we give the number of events at the LHC for an integrated luminosity of 100 
(fb)^^ of the processes metioned above as well as those of the backgound. The kinematic 
cuts listed in each column are applied sequentially. The two charged leptons in the final state 
can be either electron or muon. We take it into account by summing over the two lepton 
species. We also summed over three species of neutrinos in the ZZ background. In order 
to generalize our study, we also examine the collider reach when AniAOHO = 20 GeV and 
30 GeV. As shown in the table, a staitistical significance S/ \fB > 3 can be achieved after 
imposing the optimal cuts and it is further increased after imposing the mass window cut for 
AniAOHO = 10 GeV and 20 GeV, but slightly decreased for Attiaoho = 30 GeV because the 
mass window cut, mu < 10 GeV, also cuts away ~ 32% signal events. With an integrated 
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TABLE III: Number of signal and background events, for the benchmark points [mjjo , nij^o , mfj+] = 
[50, 60(70 , 80), 170] GeV, at the LHC with an integrated luminosity 100 fb"^ The kinematic cuts 
listed in each column are applied sequentially. Note that we summed over two species of charged 
lepton, i.e. electron and muon, for both signal and background and we also summed over three 
species of neutrinos in the ZZ background. 



BKGD 


basic 


optimal 


mu < 10 GeV 


WW 


1.1 X 10^ 


110 


62 


ZZ 


2.1 X 10^ 


3 





total 


1.3 X 10 '^ 


113 


G2 




Signal 


basic 


optimal 


mu < 10 GeV 


(50, 60) 


117 


37 


37 


S/B 


9 X lO-'' 


0.33 


0.60 


s/Vb 


0.32 


3.48 


4.70 


(50, 70) 


433 


56 


50 


S/B 


3.3 X 10-3 


0.50 


0.81 


s/Vb 


1.20 


5.27 


6.35 


(50, 80) 


680 


38 


26 


S/B 


5.2 X 10-3 


0.34 


0.42 


s/Vb 


1.89 


3.57 


3.3 



luminosity 300 ft)- , all three cases can reach 5a significance even before the mass window 
cut, therefore, LHC has a great potential to observe the dark scalar. 

We note that a competition between the basic cut and mass window cut occurs. For 
example, the acceptance of the basic cut is large but that of the mass window cut is small 
for Am^oj/o = 30 GeV. On the contrary, the acceptance of the basic cut is small but that 
of the mass window cut is large for Am^o = 10 GeV. When the mass gap is less than 
5 GeV, the pr of two charged leptons is so small that most signal events fail the basic cut, 
and hence there is no hope to observe a signal with such a small mass gap at the LHC. 
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V. CONCLUSION 



In this work we have shown how the particles of the dark scalar doublet may be discovered 
at the LHC, and how their presence will affect the properties of the SM Higgs boson. From 
present LEP data, we obtain the bound mj^^o+m^o > mz, whereas the SM Higgs mass bound 
of 114.4 GeV can be relaxed down to about 106 GeV, assuming that mA_o — m/^o = 10 GeV. 
The dark scalars dramatically affect the search of the SM Higgs boson in the intermediadte 
mass region, i.e. rrih ~ 100 — 150 GeV, at the LHC when the decay mode h H^H^ opens. 
We find that the decay branching ratios of the usual decay modes of the SM Higgs boson, 
i.e. h — s> hh/WW* / ZZ* /'-yj, are highly suppressed, ~ 60%, which makes it more challenging 
to observe the SM Higgs boson in those usual decay modes. On the contrary, one could 
detect the SM Higgs boson in its decay into the dark scalars, which will escape the collider 
detection. We show that it is very promising to search for the SM Higgs boson through its 



invisible decay in the so-called weak-boson fusion process [3J] which can cover most of the 
parameter space {rriHO, ^2) of the DSDM at the 5o" significance. After examining the decay 
pattern and decay branching ratios of the dark scalars, we consider their discovery potential 
at the LHC via the process of A^H^ associated production. For 111^0 of about 50 GeV (which 
is also consistent with its being a dark- matter candidate), it should be observable at the 
LHC. 
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